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The Fermi space telescope has detected over 100 pulsars. These discoveries have ushered in a new era of pulsar 
astrophysics at gamma-ray energies. Gamma-ray pulsars, regardless of whether they arc young, old, radio-quiet 
etc, all exhibit a seemingly unifying characteristic: a spectral energy distribution which takes the form of a 
power law with an exponential cut-off occurring between ~1 and ~10 GeV. The single known exception to 
this is the Crab pulsar, which was recently discovered to emit pulsed gamma rays at energies exceeding a few 
hundred GeV. Here we present an update on observations of the Crab pulsar above 100 GeV with VERITAS. 
We show some new results from a joint gamma-ray /radio observational campaign to search for a correlation 
between giant radio pulses and pulsed VHE emission from the Crab pulsar. We also present some preliminary 
results on Lorentz invariance violation tests performed using Fermi and VERITAS observations of the Crab 
pulsar. 



1. Introduction 

The Fermi space telescope has detected over 100 
new gamma-ray pulsars, the bulk of which are young, 
high-£/, pulsars. These discoveries have ushered in 
a new era of research into pulsars and the physics 
of emission from pulsar magnetospheres. However, 
despite the wealth of new data in recent years, the 
decades-old problem of the origin of gamma-ray emis- 
sion from pulsars remains unsolved. The origin of the 
coherent radio emission from pulsars is also poorly 
understood, as are the exotic temporal phenomena of 
Nulling and Giant Radio Pulses (GRPs) exhibited by 
some radio pulsars. 

In spite of the difficulty faced to explain the emis- 
sion from pulsars, their regular and predictable pulsa- 
tions and immense density makes them incredibly use- 
ful laboratories for the study of space-time and grav- 
itation. In this contribution, we present an update 
on observations of the Crab Pulsar above 100 GeV 
with the Very Energetic Ray Imaging Telescopes Ar- 
ray System (VERITAS). In particular, we focus on a 
recently completed study to search for a correlation 
between the emission of GRPs and pulsed very-high- 
encrgy (VHE) gamma-ray emission and preliminary 
results from a study searching for the effects of quan- 
tum gravity on the propagation of photons of different 
energies over large distances. 



2. Searching for a VHE-GRP Correlation 
2.1. GRPs in the Crab Pulsar 

The Crab pulsar, PSR B0531+21, is a powerful 
young pulsar and is, s o far, the only pulsar to be de- 
tected above 100 GeV [Aliu et al.ll201lT lAleksic et al.l 
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Figure 1: Panel (a) shows the VHE gamma-ray pulse 
profile (E-y > 1 20 GeV) of Crab pulsar measured by 
VERITAS (see I Aliu et all [20111 ). The overlaid (green) 
curve was determined from a maximum likelihood fit to 
the unbinned VERITAS data. Panel (b) shows the 
VERITAS (blue) and GRP (red) profiles for 11.6 hours 
of simultaneous VERITAS/GBT observations. 
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l2012j |. The Crab pulsar is also one of only sev- 
eral known pulsa rs exhib i ting t he phenomenon of gi- 
ant radio pulse [Knightl |2006| : single radio pulses 
with flux densities that greatly exceed the aver- 
age pulse flux density. GRPs appear in the same 
phase range as normal radio p ulses, but their en- 
ergy distribution is a po wer-law [Cordes et al.112003 , 
IPopov fe Stappersl 12007] in contrast to the regular 
pulse s that follow a Gaussian o r log-normal distribu- 
tion Burke-Spolaor et al.l 12012 1 . In the Crab pulsar, 
individual GRPs ca n be a few nanoseco nds to a few 
microseconds wide Hankins et al J [20031 ] and can, at 
their maximum, be as bright as a few million Janskys 
jSoglasnovl 20071. 

Above 6 GHz, stark differences in the behaviour of 
main pulse GRPs and intcrpulse GRPs have been ob- 
served in the Crab pulsar. Intcrpulse GRPs are typ- 
ically several microseconds long an d populate a set 
of regularly spaced frequency bands jHankins fc Eilekj 
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(a) Event Arrival Times (b) Event Separation 



Figure 2: This figure shows examples of the match between the VERITAS data and the Monte-Carlo model. The event 
rate distribution is shown on the left and the arrival time event-separation distribution for each pair of consecutive 
events is shown on the right. 



12007 in contrast to main pulse GRPs which ex- 
hibit broadband spectra and appear as a succession 
of narrow pulses ranging from unresolved widths be- 
low 0.4 ns to widths of a few microseconds. In the 
model of iLvutikovl [2007 [. the band structure seen in 
the Crab pulsar interpulsc GRPs at high frequencies 
is generated by particles in the outer pulsar magneto- 
sphere which have been energised by a magnetic field 
reconnection and emit via anomalous cyclotron res- 
onance. A prediction of this model is the existence 
of a particle beam in the outer magnetospherc with 
a large Lorentz factor, O(10 8 ), which should gener- 
ate curvature photons with energies as high as tens 
of GeV. The mechanisms responsible for the gener- 
ation of GRPs are, however, still largely unknown. 
Changes in the coherence of the plasma beam, which 
is believed to be responsible for the normal pulsed ra- 
dio emission, can in principle explain the generation 
of GRPs. Such coherence changes will, however, have 
no effect on the incoherent emission from pulsars (op- 
tical, x-ray, gamma-ray). Mechanisms which increase 
the rate of particle production within the magneto- 
spheric emission region, or which change the direc- 
tion of the emission beam can, in principle, affect the 
higher energy incoherent emission and may create an 
enhancement in the gamma-ray emission. 

Several previous studies have been performed to 
search for a possible connection between GRPs in 
the Crab and higher-energy incoherent emission. 
While only upper -limits have been report at x- 
ray (1.4-4.5 keV) jBilous et ail |2012t and gamma- 
ray (50-220 keV and 0.1-5 Ge V) iBilous et al|l201lL 
iLundgren et al.l fl995| energies. IShearer et al. [2003j . 
observed a significant 3% enhancement in the optical 
main pulse concurrent with the period of emission of 
GRPs measured at 1380 MHz. No enhancement was 
seen in the interpulse. This observation suggests a 
link between coherent radio emission and incoherent 



optical emission in the Crab pulsar. Small changes 
in the pair-creation rate leading to localised density 
increases within the emission plasma could create the 
GRP event and provide a small enhancement in the 
optical incoherent emission. 



2.2. Observations and Search Strategy 

VERITAS and the 100-m Robert C. Byrd Green 
Bank Telescope (GBT) simultaneously observed the 
Crab pulsar for a total of 11.6 hours across four 
nights in December 2008 and November and Decem- 
ber 2009. A total of 15366 GRP radio events recorded 
at 8.9 GHz, whose peak signal-to-noise ratio was at 
least seven times greater than the average radio sig- 
nal, were selected for the correlation study (see Fig- 
ure 1). The VERITAS observations were made at 
high elevation (El > 70°) under the best possible 
sky conditions with the four telescopes fully opera- 
tional and with a GPS time-stamp event tagging er- 
ror < 10 /is. The VERITAS event reconstruction and 
cosmic-ray rejec tion followed the identical procedure 
as described in [Aliu et al.l 1201 1| . The GPS time of 
the candidate gamma-ray events which passed analy- 
sis cuts was converted to barycentric dynamical time 
and phase-folde d using the Crab pulsar monthly tim- 
ing ephemeris 1 Lvne et aT1ll993l ]. 

Using the VERITAS and GRP photon lists a corre- 
lation analysis was performed. VERITAS events were 
selected if their barycentric arrival time was within 
a time window positioned around the GRP arrival 
times. Gamma-ray enhancements which lagged, lead 
or were contemporaneous with the GRP were inves- 
tigated. Since the correlation length is unknown (if 
a correlation exists at all), eight different coincidence 



1 http : //www . jb . man . ac . uk/pulsar/crab . html 
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Figure 3: 95% confidence level upper limit on the VHE 
gamma-ray flux correlated with GRPs. The slight 
negative and positive shift of the x-position of the before 
and after symbols is done as a visual aid. Limits are 
given in units of the average V HE pulsar flux measured 
by VERITAS [Aliu et al.ll20lH ]. 



limit on the enhanced VHE emission correlated with 
GRPs was calculated for each of the 72 searches per- 
formed. These upper limit values are shown in Fig- 
ure 3. 

These findings agree with res ults previously re- 
ported at keV and GeV en ergies Bilous et all l201ll 
I20121 iLundgren et al.l [l995j . GRP-emission mecha- 
nisms associated with changes in plasma coherence 
will not cause enhancements in incoherent emission, 
while small and localised changes in the pair-creation 
rate can exp lain the small ( 3%) o ptical-GRP correla- 
tion seen by IShearer et ail [2003| but yield VHE en- 
hancements below the current VHE sensitivity. In re- 
cent m odels of pulsed VHE emission from the Crab 
pulsar jAharonian et al.l I20121 iLvutikov et~aH l2012j 
any expected GRP-VHE enhancements are also be- 
low th e sensitivity of cu rrent VHE instruments. 

See lAliu et al.l [2012 [ for a detailed description of 
this VHE-GRP correlation study. 



3. Searching for Quantum Gravity 



time windows were chosen lasting; 1, 3, 9, 27, 81, 243, 
729 and 2187 pulsar rotations. Main pulse and in- 
terpulse GRPs were also considered separately, thus 
72 different correlation searches were performed. A 
Monte-Carlo model of the VERITAS data was used 
to determine the number of coincident VHE events 
expected in the absence of a correlation. Monte-Carlo 
time series data sets were generated by drawing ran- 
dom times from the raw VERITAS trigger rate his- 
tograms. A random subset of these events was chosen 
to contain the temporal signature of the pulsar with 
phase values draw from the functional form which de- 
scribes the shape of the pulse profile of the Crab pul- 
sar as measured by VERITAS (see the dashed curve 
in panel (a) in Figure 1). This procedure was used 
to generate sets of Monte-Carlo data which contain 
a signal from the Crab pulsar at a chosen flux level 
and which model all the data rate characteristics of 
the real VERITAS data. Examples of the match be- 
tween the Monte-Carlo time-series datasets and the 
VERITAS gamma-ray data are shown in Figure 2. 



Many theoretical investigations aim to formulate a 
theory of quantum gravity in order to unify the the- 
ories of quantum mechanics and general relativity. 
Some approaches to the theory of quantum gravity 
result in equations of motion which contain an energy 
dependent dispersion relation for massless particles, 
meaning that the speed of light in a vacuum would 
not be constant. Such theories would violate Lorentz 
invariancc. Depending on the particular framework, 
photons would travel a the speed 



v(E) = c 



1 ± 



E 



E, 



QG 



(f) 



where E is the photon energy and Eqq is the energy 
scale at wh ich the effects of quantum gr avity become 
important Amelino-Camelia et al.lll998| . The order, 
n, of the energy dependence, and the ± sign ambiguity 
would be fixed within the framework of a given model 
for quantum gravity. To search for these dispersive 
quantum gravity effects one can apply the following 
framework: 



2.3. GRP-VHE Correlation Results 

No significant enhancement in VHE gamma-ray 
emission correlated with GRPs measured at 8.9 GHz 
was observed in any of the 72 searches performed. Us- 
ing the Monte-Carlo model of the Crab pulsar signal 
within the VERITAS data, it is possible to investi- 
gate the probability that a given number of correlated 
events is consistent with a chosen level of correlated 
VHE enhancement . Following a similar procedure as 
iBilous et al. I [201l| . the 95% confidence level upper- 



• Two photons with energies E\ and Ei emitted 
simultaneously arrive at an observer at a dis- 
tance, L, with time separation, At. 

• Depending on the order n of the energy depen- 
dence, the time difference At for the linear and 
the quadratic terms in {E/Eqq) are: 



. L Eo — E\ 
Ah = - 

c E QG 



E QG = 2 — -± 2 

c A<i 



and 
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L3EZ-E? „ IL3E2-E? , n , 

• By measuring the structure in a temporal event, 
such as an emission flare, in different energy 
bands, we can place an upper limit on At. 

• By placing a limit on the maximal At, a lower 
bound on Eqg can be found. 

From Equations 2 and 3 it is clear that when search- 
ing for Lorcntz invariance violation (LIV) it is favor- 
able to observe temporal phenomena which occur on 
short time scales (minimising At) at a large distance 
from the observer (maximising L) over a wide pho- 
ton energy range (maximising E2 — E\). Gamma-ray 
bursts (GRBs) and TeV gamma-ray flares from active 
galactic nuclei (AGN) meet all of the above criteria, 
and to date, have been used to place the strongest lim- 
its on the energy scale of quantum gravity [H.E.S.S. 
Collaboration ct al. 2011, Abdo ct al. 2009]. 



3.2. Preliminary Results 

Using the VERITAS and Fermi Crab pulse profiles 
shown in Figure 4 we can place a limit on Eqq. The 
positions of the emission peaks are consistent within 
measurement uncertainty. The measurement of the 
VERITAS emission peak positions, which was deter- 
mined from an unbinned maximum likelihood fit to 
the VERITAS phase data (solid line plotted in Fig- 
ure 4), is the dominant source of uncertainty. This 
uncertainty is about 2 x 10~ 3 phase units, or, 60 /is. 
Using the 95% confidence level upper limit on the time 
separation between the peaks, 100 (is, the canonical 
distance to the Crab pulsar of 2 kpc, and the con- 
servative value of 120 GeV for the difference between 
the Fermi and VERITAS energy ranges, we arrive at 
3 x 10 17 GeV and 7 x 10 9 GeV for the lower limit on 
the values of Eqg for the liner and quadratic case, 
respectively. 



4. Conclusion 



3.1. LIV with Gamma-ray Pulsars 



iKaaretl fT999j showed that gamma-ray pulsars were 
also excellent tools for LIV searches and achieved a 
lower limit on the energy scale of quantum gravity 
of 1.8 x 10 15 GeV for case of re = 1 (Equation 2). 
Further, doing LIV tests with gamma-ray pulsars has 
several benefits over flaring transients such as GRBs 
or AGNs: 

• Tests do not rely on random transient events 
and observational luck. 

• Limits improve with longer exposures. 

• Limits are based on highly significant measures 
of peak positions rather than on a handful of 
single photons. 

• Delay effects intrinsic to the source can be dis- 
tinguished from LIV effects. 

• If a delay effect is measured, the delay separa- 
tion will appear as a constant phase offset as the 
pulsar period lengthens, if it's a source effect. 

• If a delay effect is measured and is caused by 
LIV, this effect should stay constant in time as 
the pulsar slows and so should shift in phase by 
the amount A$(i) = At/(P + tP). 

Following the recent detection of the Crab pulsar by 
VERITAS, which has extended the observed energy 
range of the pulsar by over an order of magnitude, 
we are promp ted to repeat the LIV test performed by 
IKaaretl j!999j . 



VERITAS is engaged in a broad and varied pul- 
sar observation campaign. These proceedings present 
a summary of two research lines on the Crab pulsar 
which have recently been, or are soon to be, com- 
pleted. These include the first limits on the level of 
enhanced VHE emission correlated with GRPs in the 
Crab pulsar. Given the level of uncertainty in the- 
ories of GRP emission, it is hard to draw firm con- 
clusions resulting from the lack of any observed cor- 
relation with VHE emission. We are not aware of 
any theory with quantitative predictions of correlated 
emission between GRPs and VHE emission. Small 
and localised changes in the pair-creation rate, which 
can explain the smal l (3%) optical enhan cements pre- 
viously measured bv lShearer et al.l |2003l |. would yield 
VHE flux enhancements which are below our sensitiv- 
ity. 

We also present some preliminary limits on the en- 
ergy scale of quantum gravity determined from VERI- 
TAS and Fermi observation of the Crab pulsar. These 
limits are only one order of magnitude lower than the 
best limits ever achieved (with GRBs and AGNs) and 
are set to improve with more data, more VHE pulsars, 
and improved analysis methods. 
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(a) Crab pulsar main pulse (b) Crab pulsar intcrpulse 

Figure 4: A zoom in of the Crab pulsar gamma-ray emission peaks [Aliu et al ] |201lj . The full VERITAS profiles are 
shown in Figure 1. T he overlaid curve i s from a maximum likelihood fit to the unbinned VERITAS data. The Fermi 
points are taken from lAbdo et al] [20101 ]. 
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